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A first-principles approach to the construction of concentration-temperature magnetic phase diagrams of 
metallic alloys is presented. The method employs self-consistent total energy calculations based on the co¬ 
herent potential approximation for partially ordered and noncollinear magnetic states and is able to account for 
competing interactions and multiple magnetic phases. Application to the Fei_^MnjPt “magnetic chameleon” 
system yields the sequence of magnetic phases at T = 0 and the c-T magnetic phase diagram in good agree¬ 
ment with experiment, and a new low-temperature phase is predicted at the Mn-rich end. The importance of 
non-Heisenberg interactions for the description of the magnetic phase diagram is demonstrated. 


Magnetic substitutional alloys are often found to excel in 
applications mil, because alloying broadens the parame¬ 
ter space for tuning the desired properties. However, wide 
tunability, combined with the need to target certain operating 
temperature ranges, presents a challenge for empirical mate¬ 
rials design. Competing magnetic interactions in alloys can 
produce complicated magnetic phase diagrams (MPD) with 
multiple magnetic phases ffl. Understanding of the c-T 
MPD’s is therefore essential for the development of advanced 
magnetic materials. 

Some MPD’s can be computed using the Heisenberg model 
with empirical or calculated exchange parameters combined 
with the mean-field approximation (MFA) ItI I^ M onte-Carlo 
simulations iMl, or spin-fluctuation theory ifoll . However, 
many systems are not adequately described by the Heisenberg 
model. In metallic alloys the interaction parameters are sensi¬ 
tive to the electronic structure and population and thereby to 
the content of the alloy ||il0 n ] and to the degree of spin dis¬ 
order El. To avoid the limitations of the Heisenberg model, 
one can use first-principles spin-dynamics simulations 01511 
or construct a generalized spin Hamiltonian to map the adi¬ 


abatic energy surface 0161 llTll for use in thermodynamic cal¬ 


culations. The energies of disordered spin configurations can 
also be obtained using the disordered local moment (DLM) 
model |T§, 19], where the spin-rotational averaging is done in 
the coherent potential approximation (CPA). While all these 
approaches fail in strongly itinerant magnets, they are applica¬ 
ble when the spin moments do not vary by more than 10-20% 
in different spin configurations. We restrict ourselves to such 
systems here. 

First-principles spin-dynamics and the construction of a mi¬ 
croscopic generalized spin Hamiltonian are computationally 
very demanding and unfeasible for most systems of practi¬ 
cal interest. We have developed an alternative approach, 
in which self-consistent DLM and noncollinear CPA calcu¬ 
lations are used to construct a Ginzburg-Landau-type total- 
energy functional expressed through a small number of rele¬ 
vant magnetic order parameters. Combined with the MFA ex¬ 
pression for the magnetic entropy, this method provides the 
variational free energy. A similar scheme was used to de¬ 


scribe the phase transitions in FeRh 0211] . Here we show, 
using the Fei-^Mn^^Pt disordered alloy system as a test case, 
that this efficient approach is sufficiently powerful to explain 
and refine a complicated MPD, not only reproducing the five 
known magnetic phases but also predicting another, hidden 
low-temperature phase in this system. 

Fei^jMn;cPt alloys are of interest for ultrahigh-density mag¬ 
netic recording and medical applications . Their structural 
ordering is of the LIq type in the fee sublattice, with (001) 
layers of Pt alternating with disordered Fe/Mn layers. Neu¬ 
tron diffraction measurements revealed three collinear and two 


noncollinear phases Il23n . The collinear phases are the ferro¬ 
magnetic (F) phase at the Fe-rich end, the C-type antiferro¬ 
magnetic phase at the Mn-rich end, and the G-type antiferro¬ 
magnetic phase in the middle of the diagram. The correspond¬ 
ing ordering wave vectors are Qf - (0,0,0), Qc = (1,0,0), 
and Qg = (1,0,1/2) in units of Inja (or Injc for the z com¬ 
ponent). The transitions between the collinear phases occur 
through intermediate 2Q phases combining the corresponding 
orderings for two orthogonal spin components (see Supple¬ 
mental Material 0 for an illustration). 

The samples studied in Ref. l23| show a high degree of LIq 
order at all concentrations, which is consistent with the fact 
that magnetic ordering occurs well below the structural or¬ 
dering temperatures in these alloys. Therefore, as a practical 
simplification we assumed perfect LIq ordering and complete 
disorder of Fe and Mn atoms within their own sublattice. The 
method can also be applied to alloys with partial chemical or¬ 
dering and in principle allows one to study the coupling be¬ 
tween magnetic and chemical order parameters 11251] . 

To construct the energy functional, we have extended our 
CPA implementation ||26 


27l] based on the tight-binding lin¬ 


ear muffin-tin orbital formalism Q by special features de¬ 
signed to describe complicated magnetic states. First, we im¬ 
plemented the vector DLM (VDLM) model, in which partially 
ordered magnetic states are specified by the Curie-Weiss dis¬ 
tribution functions pifi{d) oc exp(a/^ cos 6), where i is the lat¬ 
tice site and p the component index; are regarded as vari¬ 
ational parameters. This formulation is suitable for systems 
with axial spin symmetry, i. e. those with collinear magnetic 
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order. Different spin moment orientations are treated as dif¬ 
ferent alloy components in the CPA formalism. The integral 
over the azimuthal angle in the CPA equations is taken analyt¬ 
ically, while the 6 dependence is discretized using the 16-point 
Gauss-Legendre quadrature. The potentials for all atoms are 
determined by embedding the CPA self-consistency loop into 
the density-functional theory (DFT) charge iteration. To en¬ 
force magnetic self-consistency, constraining transverse mag¬ 
netic fields Ellii are introduced for each orientation of the 
spin moment and determined self-consistently. Local density 
approximation is used for exchange and correlation. 

The second feature extends CPA to the noncollinear case, in 
which the orientations of the spin moments of different com¬ 
ponents on the same lattice site can be different. This method 
is suitable, in particular, for the studies of 2Q structures ap¬ 
pearing in the Fei_j:MnxPt system. Self-consistent constrain¬ 
ing fields are also used in these calculations. Both VDLM and 
noncollinear CPA calculations yield the DFT total energy. 

Let us first examine the magnetic interactions in the vicin¬ 
ity of the paramagnetic state. We set up a unit cell for each of 
the three magnetic orderings and calculate the total energy for 
about 70 partially ordered VDLM states with lapel and laMnI 
ranging from 0 to 3. Experimental lattice constants are used 
at the concentrations reported in Ref. At each concen¬ 
tration the calculated total magnetic energy Emag (per Fe/Mn 
atom, referenced from the paramagnetic state) is fitted to 
even fourth-order polynomials in the reduced magnetizations 
= {cos 6ft), which are expressed through the fields ay by 
the Langevin function. The quadratic part of these polynomi¬ 
als, Eq = 57FeFe(Q)'«Fe+5-^MnMn(Q)»JMn + -^FeMn(Q)»JFe»JMn, 
defines the component-resolved effective exchange interac¬ 
tions JftviQ) for the three orderings. Since the total energies 
are calculated with constraining fields, these exchange param¬ 
eters are free from the errors associated with the long-wave 
approximation Bill . The results are shown in Fig.[T] 

The concentration dependence of the parameters 7^y(Q) 
shows that, in agreement with experiment, the F, G, and 
C-type orderings are favored at the Fe-rich end, in the 
middle, and at the Mn-rich end, respectively. Further in¬ 
sight can be obtained from the reduced exchange param¬ 
eters that are normalized by the concentrations, JftviQ) - 
JiiviQ)liCftCv), shown in Fig. |2l These reduced quantities 
would be concentration-independent in a Fleisenberg system 
with pair exchange parameters 7^v(R) depending only on the 
distance and the identity of the atoms in a pair. We see that the 
parameters JftviQ) for like components (i. e. Fe-Fe and Mn- 
Mn) are almost constant for all ordering vectors, as well as 
the reduced Fe-Mn coupling for the F ordering. However, the 
reduced Fe-Mn couplings at the G-type and C-type ordering 
vectors depend strongly on the concentration, which reflects 
the effect of band filling on the exchange interaction in metal¬ 
lic systems. 

The local spin moments of Fe and Mn increase by about 
10% as X goes from 0 to 1 (for example, from 2.87 to 3.11 
for Fe and from 3.35 to 3.74 for Mn in the paramagnetic 
phase), and similar variations are observed for different phases 


and spin directions. The adiabatic approach S is thus well 
suited for this system. We also repeated some calculations 
without the constraining magnetic fields, which is equivalent 
to making the long-wave approximation iImIi . and found that 
the resulting errors in T^y for all phases do not exceed 5-7%. 

To further extend the mapping of the magnetic configura¬ 
tion space, we performed noncollinear CPA calculations for 
the F/G, G/C, and F/C noncollinear 2Q phases in the relevant 
concentration ranges. These calculations are needed to reveal 
the possible interaction between orderings at different Q in the 
2Q structures, which can appear in quartic and higher-order 
interaction terms. A 2Q structure is parameterized by two an¬ 
gles, 0Fe and Amii that the spin moments of Fe and Mn atoms 
make with the z axis. (The z and x components of the magne¬ 
tization order with one or the other of the Q vectors.) Using 
the symmetries, the accessible space of (0Fe, ^Mn) is reduced 
to the range 0 < 0Fe ^ ^/2, -n < 0Mn ^ with additional 
^Mn —>?!■- ^Mn Symmetry at 0Fe = ■7!'/2 and Oyin -^Mn 
symmetry at 0Fe = 0. This irreducible domain is covered by a 
uniform mesh of 38 points. 

For a particular 2Q phase (say, F/G) we then combine the 
data from the separate VDLM calculations for the F and G 
phases with those from the noncollinear CPA results for the 
2Q phase and fit the magnetic energy at the given concen¬ 
tration to a polynomial in the order parameters mFe.F, '«Fe,Gj 
wiMn.F, and mjvin.G- We allowed all symmetry-respecting terms 
in the poly nomials of up to sixth order (see Supplemental Ma¬ 
terial 1241] for details). The high accuracy of the fits is illus¬ 
trated in Fig. for the F/G phase at .r = 0.26; all other fits 
are of similar accuracy. Fig. [Jr shows the magnetic energy 
as a function of 0Fe, ^Mn in the same phase at x = 0.26. At 
this concentration there are two minima: the global one at 
(7r/2,7r/2) corresponding to the collinear G-type phase, and a 
local one near in 16, 7r/2) corresponding to the 2Q phase. At a 
lower concentration there is a first-order transition where the 
2Q minimum goes below the G-type minimum. 

Using the combined fits for the magnetic energy, we now 
determine the ground states for all concentrations by choosing 
the lowest energy of all the competing phases. The results are 
shown in Fig.|4] We find that the F phase is stable at x < 0.11 
and the G phase at 0.23 < x < 0.66. There is a 2Q F/G phase 
at 0.11 < X < 0.23 separated from F by a second-order and 
from G by a first-order transition. The G/C phase is stable at 
0.66 < X < 0.85. 

Surprisingly, at x = 0.85 we find a first-order transition 
from the G/C to the F/C phase. The existence of this first-order 
transition is in excellent agreement with the observed a brup t 
drop in the mean magnetic moment at this concentration 12311 . 
However, since the F component in the low-temperature Mn- 
rich phase has not been previously identified, the existence of 
the F/C phase is a prediction that needs to be verified experi¬ 
mentally. The F/C and G/C phases at the Mn-rich end differ 
essentially in the Fe ordering alone, as the ordering of the Mn 
spins is almost purely C-type. The energy difference between 
the F/C and G/C phases reaches about 20 meV per Fe atom 
near x = 0.95; it is barely visible in Fig. |4] because of the 
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FIG. 1. Magnetic exchange parameters in the paramagnetic state corresponding to (F) ferromagnetic ordering, Qf = 0; (G) G-type 
ordering, Qg = (1,0,1/2); (C) C-type ordering, Qc = (1,0,0). Red circles: ypePeCQ); blue diamonds: yMnMn(Q); black squares: ypcMn(Q). 
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FIG. 2. Same as in Fig.[T]but normalized by the concentrations: y^v(Q) = ■4iv(Q)/(Cf,Cv). 
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FIG. 3. (a) Accuracy of the fit at x = 0.26. (b) The magnetic energy 
Emag (meV per Fe/Mn atom) in the F/G phase at x = 0.26. The global 
minimum at {n/l, n/2) corresponds to the collinear G-type phase. 



small Fe concentration. 

The first-order transition from F/G to G at x » 0.23 is also 
in excellent agreement with experiment, while the transition 
from G to G/C occurs at a larger x compared to experiment, 
where it is close to x = 0.5. Note, however, that our calcu¬ 
lations are for systems with perfect 3dl5d ordering, while the 
order parameter in the experimental samples for x = 0.5 and 
x = 0.6 was 0.79 Q. 

Fig. |5] shows the angles Ope, Amh in the ground states as a 
function of x. (By convention, F and C amplitudes lie along 
the z and G along the x axis; except for the F/C phase, where 
F is along x.) As could be inferred from Fig. [T] the spin 
moments of Fe and Mn are antiparallel in the F phase and 
parallel in the G phase. First-order phase transitions appear 
as discontinuous jumps of the angles. Note that in a wide 


FIG. 4. Magnetic energies of different phases at zero temperature 
(per Fe/Mn atom). Arrows show the boundaries between the F, F/G, 
G, G/C, F/C, and C phases (in the order of increasing x). 


concentration range at the Mn-rich end the Fe spin moments 
are almost perpendicular to the MnPt host in the ground state, 
while Mn ordering is almost pure C-type. (A full set of first- 
principles calculations with interpolated lattice constants was 
performed with a small 0.01 step in 0.94 < x < 1 range to 
confirm this.) This feature highlights the importance of non- 
Heisenberg terms in the magnetic energy. 

We now turn to the full concentration-temperature phase 
diagram. We have the fits for the magnetic energy 
Fmag(niFe, DiMn), where the x and z components of mpe and 
iiiMn correspond to the two Q vectors; in the collinear phase 
one of these components vanishes. Emag is available for dis- 
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FIG. 5. Angles By made by the spin moments of Fe and Mn with the 
z axis at zero temperature. By convention, the G order parameter, as 
well as the F order parameter in the F/C phase, are assumed to be 
orthogonal to the z axis; all others are parallel to z. 


Crete concentrations, and we use linear interpolation between 
them. The entropy is approximated as S = (1 - x)S(mps) + 
xS {niMn), where S (m) is the entropy of a classical spin in a 
Weiss held of such a magnitude that the magnetization is m. 
This corresponds to the mean-held-like distribution function 
Py{ff) oc exp(av cos ff), where ff is the angle with respect to 
the direction of niy. Given that mpe.F = '«Fe cos 0Fe in the F/G 
phase, etc., we minimize the free energy for each phase with 
respect to four parameters 0Fe, ^Mn, o^Fe and aMn- 


For a system with purely Heisenberg interaction this 
scheme is identical to MFA applied to 2Q phases fl. There¬ 
fore, both second-order transitions to the paramagnetic phase 
(where non-Heisenberg terms have no effect in MFA) and 
ground-state properties (where only energy is important) are 
correctly described by this approach. In the intermediate 
temperature range our scheme can therefore be treated as a 
thermodynamic interpolation. At such intermediate temper¬ 
atures we are essentially assuming that the presence of non- 
Heisenberg terms does not strongly change the distribution 
functions and that our magnetic energy fits remain valid for 
partially disordered 2Q phases. In most cases, these assump¬ 
tions are likely consistent with the accuracy of the MFA it¬ 
self. Significant improvement of the thermodynamic descrip¬ 
tion could be reached by mapping the total energies to a mi¬ 
croscopic spin Hamiltonian, but this expensive procedure is 
beyond the scope of this paper. 


Fig.|6]shows the MPD obtained both from the full magnetic 
energy hts and from the same hts truncated at the quadratic 
(Heisenberg) terms. First-order transitions are shown by 
dashed lines. The overall structure of the phase diagram 
agrees well with experiment. The first-order transition from 
the F/G to the G phase may help explain the peculiarities 
of the observed temperature-dependent magnetization peaks 
I 23n . Indeed, there is a concentration range from x - 0.23 to 
X X 0.28 where the ground-state G-type ordering turns into 
F/G and then to F on heating. The physics is complicated by 


conhgurational disorder, which may lead to the formation of 
small Fe-rich clusters at elevated temperatures 11231] . 

If only Heisenberg terms are kept in the magnetic energy, 
the transitions into the paramagnetic phase remain unchanged. 
However, all first-order transitions turn into second-order; the 
F/G to G transition is shifted to much larger concentrations, 
and the F/C phase disappears completely. Thus, in order to 
describe the observed first-order transitions at x - 0.23 and 
X = 0.85 it is important to take the non-Heisenberg interac¬ 
tion terms into account. The F/C phase disorders at rather low 
temperatures in our description, although its stability may be 
underestimated in our thermodynamic scheme. The non-zero 
magnetization in this phase should facilitate an easy experi¬ 
mental verification of its existence. 



FIG. 6. Magnetic phase diagram of Fei_jMnj^Pt. Temperatures are 
rescaled by the ratio Tgyp(x)/T,i,(x), where T(x) = (1 - x)Tc + xT^; 
Tc and Fa; are the ordering temperatures of FePt and MnPt from Ref. 
1^ {exp) or theory (ih). (MFA gives Tc = 924 K and Tpi = 1670 
K.) Thick (blue) lines: full fit for Thin (red) lines: same fit but 
with non-Heisenberg terms set to zero. Solid (dashed) lines: second 
(first) order phase transitions. Symbols: experimental data 


To conclude, we developed a computational tool based on 
a combination of first-principles calculations that is capable 
of describing a complicated c-T magnetic phase diagram of 
a metallic alloy with competing interactions. Its application 
to the Fei_j:MniPt system produced a detailed interpretation 
of the experimental phase diagram and also predicted the ex¬ 
istence of a previously unknown low-temperature magnetic 
phase on the Mn-rich end. The correct first-order transitions 
and the Mn-rich phase are only captured if non-Heisenberg 
terms are included in the magnetic energy, showing the limi¬ 
tations of the conventional approach based on the Heisenberg- 
model Hamiltonian. The wide applicability and predictive 
power of this approach makes it useful for the design of mag¬ 
netic materials with desired properties. 
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SUPPLEMENTAL MATERIAL 


In its basic features, our implementation IS 111 of the co¬ 
herent potential approximation (CPA) follows the formulation 
I S2t] based on iterating the “coherent interactor” matrix Q for 
sites with chemical and/or spin disorder, which is defined as 


n(z) = p(z) - g-\z), 


( 1 ) 


where z is the complex energy, P is the coherent potential and 
g the on-site auxiliary Green’s function. The latter is com¬ 
puted as the on-site block of the k-integrated crystal Green’s 
function g(k) = [Pb^bb' - S BB'(k)]“^ where B and B' denote 
basis sites within the unit cell, and S is the structure constant 
matrix of the tight-binding linear muffin-tin orbital (LMTO) 
method. The coherent potential is in turn calculated from the 
LMTO potential parameters Pa and the Q. matrix: 


P^[j]pa(Pa-nr^ 

a 


-1 


( 2 ) 


where pa is the probability to hnd an atom of the component 
a on the given site. 

In the vector disordered local moment method (VDLM), 
different directions of the local magnetic moments are treated 
as different alloy components. The component index a in Eq. 
(|2|i then combines the index of the chemical component p and 
the orientation of the local moment h^, and the probabilities 
become pa - c^Pij{h^), where is the concentration of the 
component p on the given site, and Pij{h^) is the distribution 
function for the spin orientations. 

The LMTO potential parameter matrix for each spin direc¬ 
tion is initially computed in the reference frame in which the 
polar axis is aligned with the direction b of the total effective 
magnetic held for this direction (more on this below). Let us 
denoting this matrix Pa- To obtain Pa in Eq. (|2|i, we need to 
rotate Pa to the global reference frame: 

Pa = U{b)PaU-\b) (3) 


where U{b) is the rotation operator that rotates the z vector to 
b. In this work spin-orbit coupling is neglected. Therefore, 
Pa is spin-diagonal and the rotation U(b) only affects the spin 
indices. If spin-orbit coupling is included ISJ], U{b) acts on 
both spin and orbital indices and is generated by the total an¬ 
gular momentum operator J. 

The summation over a in Eq. (|2|i involves a sum over 
chemical components and an integration over the sphere for 
each component treated within the VDLM scheme. Eor our 
present purposes the VDLM scheme is implemented for par¬ 
tially ordered magnetic states with a global axial symmetry 
and without spin-orbit coupling. This symmetry implies that 
the LMTO potential matrices Pa and the spin distribution 
functions p^{h^) depend only on the chemical identity p and 
the polar angle 6^ with respect to the magnetic symmetry axis. 
In this special case the integral over the azimuthal angle in Eq. 
|2]has a very simple effect: it eliminates all spin off-diagonal 


elements and leaves spin-diagonal elements (which do not de¬ 
pend on 0) unchanged. The remaining integration over the po¬ 
lar angle 9^ is discretized using the 16-point Gauss-Legendre 
quadrature. Thus, formally we are dealing with 32 CPA com¬ 
ponents on the Fe/Mn sites corresponding to 16 values of the 
polar angle for Fe and Mn. 

The spin distribution functions for partially ordered VDLM 
states are assumed to have the Curie-Weiss form p^(fl) oc 
exp(a^ cos 9). The factors are treated as variational param¬ 
eters in the free energy minimization. 

In addition to VDLM states, in the construction of the 
energy functional we also use CPA calculations for the 2Q 
phases (F/G, G/C, and F/C) with dehnite (non-random) spin 
directions. These calculations are essential to capture the 
higher-order interactions between the order parameters cor¬ 
responding to different Q vectors; indeed, axially-symmetric 
VDLM calculations represent partially ordered states with 
only one ordering Q vector. In these calculations, the sum 
over a in Eq. (|2]i has only two terms for Fe and Mn with their 
spins pointing in certain directions within the xOz plane (as is 
sufficient for the 2Q states). Of course, the spin off-diagonal 
terms are retained in Eq. whereby P, fi, and g also be¬ 
come spin off-diagonal. The magnetic configuration in a 2Q 
phase can be written as exp(/QiR/) -i-m^ exp(/Q 2 R/) 

with orthogonal order parameters (m^nij = 0), where R, is 
the coordinate of site i. Because of the orthogonality, there 
are four real order parameters in a 2Q phase, which can be 
simply denoted as and m^. represents the averaged lo¬ 
cal moment of an atom of component p appearing at lattice 
site i. The site magnetizations are defined as - (cos 6(,), 
i. e. the dependence of the local moment on its orientation is 
neglected in this definition. This choice does impair the accu¬ 
racy of the calculations, because the htting of the calculated 
total energies and the calculation of the magnetic entropy are 
consistent with this specihc dehnition of the order parame¬ 
ters. The relevant ordering wave vectors are = (0,0,0), 
Qc = (1,0,0), and Qc = (1,0,1 /2) in units of2nla (or 2nlc 
for the z component). All these wave vectors are commensu¬ 
rate with a tetragonal a x ax2c unit cell, which is therefore 
used for the calculation of the total energies in all magnetic 
phases. (The use of the same unit cell improves the accuracy 
of the calculations.) The magnetic structures of the consid¬ 
ered 2Q states are shown in Fig.[ST]along with the G-type and 
C-type orderings (the ferromagnetic state is not shown). 

Using the symmetries 9^^ ^ n - 9f^ and ^ (applied 
to all p at the same time), the accessible space of (0Fe,6*Mn) 
is reduced to the range 0 < 0Fe ^ zt/ 2, -n < 0Mn ^ with 
additional 0Mn —> - 0Mn symmetry at 0Fe = :^/2 and Amii ^ 

-^Mn symmetry at 0Fe = 0. This irreducible domain is covered 
by a uniform square mesh of 38 inequivalent points. 

The LMTO charges are computed in the usual way from the 
conditionally-averaged physical Green’s function iS2|] Ga - 

+ PagaPa in the reference frame in which the polar axis is 
aligned with the prescribed spin direction h. Here Aa and pa 
are the spin-diagonal potential parameters related to Pa, and 
ga - [Pa - U^^(h)Q.U(h)]^^ . We emphasize that the LMTO 
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FIG. SI. Magnetic structures: (a) G-type, (b) C-type, (c) F/G, (d) G/C, (e) F/C. Large blue (small gray) spheres show Fe/Mn (Pt) sites. Green 
and red arrows in panels (c)-(d) show the spins of Fe and Mn atoms, which have similar ordering patterns but are not parallel to each other. 
Since spin-orbit coupling is disregarded, only relative angles between the spins are important. 


charges and potentials are different for different orientations of 
the local moment, and self-consistency is achieved indepen¬ 
dently for all of them. Since the coherent interactor Q, which 
includes the magnetic coupling to the rest of the crystal, and 
Pa, describing the local exchange-correlation field, are diag¬ 
onal in different reference frames, the magnetic moment ob¬ 
tained from Ga deviates from the direction b of the total effec¬ 
tive magnetic field. If the transverse component of this output 
magnetic moment is discarded, the obtained magnetic state is 
not self-consistent in the sense of the density functional theory 
(DFT). In order to remedy this problem, constraining trans¬ 
verse magnetic fields need to be introduced llS4l ISSll. We im¬ 
plement these fields in the form suggested in Ref. |Sa Specif¬ 
ically, the constraining held for each atom is assumed to have 
the same radial dependence as the exchange-correlation held 
on the same atom. The total effective magnetic held is then 
collinear within the atomic sphere, oriented at an angle 6 to the 
exchange-correlation held, and rescaled by a factor Vl + 
iH. At each iteration toward self-consistency, the misalign¬ 
ment angle 6 is adjusted to eliminate the deviation of the mag¬ 
netic moment from the prescribed direction h. Thus, at self- 
consistency, the local moment is collinear with the exchange- 
correlation held direction h, while the total effective magnetic 
held is rotated by an angle d with respect to h. 


For the 2Q states there is an additional complication that the 
orientations of the induced local moments on the Pt atoms are 
not known in advance. These orientations are also updated 
at each iteration using the output directions of the induced 
magnetic moments, guaranteeing full DFT self-consistency in 
the hnal state. 

As explained in the main text, the VDLM and noncollinear 
CPA calculations were performed for a number of concen¬ 
trations with experimental lattice parameters. The calculated 
(about 200) total energies for each candidate 2Q phase at 
the given concentration are then htted to a polynomial in the 
four order parameters m^, (where 1 and 2 label the or¬ 
dering types, e.g. F and G for phase F/G): ^(mpmj) = 
Pi(ot^) + P 2 iin^) + P\ 2 {m^^,nf^). The polynomials Pi and P 2 
include all even terms up to sixth order, and P 12 includes all 
products of second-order monomials in and (lower- 
order terms in P 12 are forbidden by symmetry). This fit also 
covered the two 1 Q phases, which are obtained at =0 or 
= 0 . 

Fig. 3 of the main text illustrates the quality of the fit and 
its energy profile for fully ordered F/G phase at one selected 
concentration x = 0.26. The actual expression is^ is the 
following {xj = m^, yj = m“"): 

















































PF(xF,yF) — - 64.78x^ + WSAxFyF + 11.96y^ 

- 11.79x^ + Q.%5xlyF + 5.66xj-yj- + 1.79xFyF + 

+ 14.39X® + l.lSx^yF - 0.69xjyj- - 2.25x^y^ - 5.48x^y^ - 2.08xFyF " 2.00y^ 
Pc(xc,yc) = - 25.25x^ - AS.lSxcyc - 9-18^1 

+ 13.02x^ + 2.\6xlyc - - 3.16x0^0 + l-14yG 

- 17.68xS - 5.56x^yG + 0.82x^y| - 2.16xlyl + 2.224yG + 3 . 65 xg 4 " 1-24 
PFG(xF,yF, Xc,yc) = + 23.09x^Xg + 0.92xj-xcyG - 1-82x^4 + 5.72xFyFX^c + 5.42xFyFXGyG - A.SxFyry^ 

+ 9.81xg 4 + l-6l4-^GyG - 8.6144 


The polynomial P 12 describes interaction between the order 
parameters at two different Q vectors. Information needed to 
fit this polynomial is only contained in the noncollinear CPA 
results for the given 2Q phase. (P 12 vanishes in the 1 Q phases 
described by the VDLM method.) Since these calculations 
(contrary to VDLM ones) are done only for magnetically- 
ordered phases, certain higher-order terms in P 12 would be 
linearly dependent within the available dataset and could not 
be distinguished. Indeed, in a fully-ordered 2Q phase we have 
+ (m^)^ = 1, and terms like x^x^, x^x^, and XjX^ 
linearly dependent. Therefore, for most concentrations we re¬ 
stricted P 12 to fourth-order terms, which are expected to be 
the most important. However, at the Mn-rich end the interac¬ 
tion of the Fe spins with the Mn host has a particularly strong 
non-Heisenberg character in the F/C phase, making the fourth- 
order fit for P 12 insufficient. Therefore, at x > 0.8 we added 
sixth-order terms to Ppc that are second-order in F and fourth- 
order in C order parameters (such as x^x'^, etc.). In addition, 
to provide a consistent fit for both G/C and F/C phases, in 
this concentration range we combined all input data to ob¬ 
tain a single fit for the 3Q structure involving all three order 
parameters. (Possible interaction between the F and G order 


parameters, which is irrelevant in this region, was left unde¬ 
termined.) 
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